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ABSTRACT 


A pert\irbed orbit, resonant scattering theory for pitch-angle 
diffusion in magnetostatic turbulence Is slightly 

generalized and then utilized to compute the diffusion coefficient for 
spatial propagation parallel to the mean magnetic field, . All 
divergences inherent in the quasilinear formalism when the power 
spectrum of the fluctuation field falls off as k ^ (q s 2) are removed. 
Various methods of computing are ct^npared and limits on the validity 
of the theory eire discussed. For 1 < q < 2 the various methods give 
roughly conq)arable values of , but use of perturbed orbits systematically 
results in a somewhat smaller K than one obtains from quasilinear theory. 
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I. INTRODUCTION 







The propagation of charged particles throijgh interstellar and 
interplanetary space has often been described as a- random process in which 
the particles are scattered by ambient electromagnetic turbulence. In 
general, this changes both the magnitude and direction of the particles’ 
momentum. In this T^aper we confine our attention to those situations 
for which scattering in direction (pitch angle) is of primary interest. 

One can derive from this microscale phenomenon of pitch-angle scattering, 
a spatial diffusion coefficient, , which describes the macroscopic 
diffusion of the particle distribution. These diffusion coefficients 
liave been useful for describing the solar modulation of the galactic 
cosmic radiation, in which the outward convection due to the solar wind 
is balanced by an inward diffusion [v. Jokipii . 1971 for a review]. In 
addition, particle observations during solar flares often show a diffusive 
phase in which the initial anisotropy of the particle distribution decays to 
isotropy. Again, knowledge of aids in calculating various characteristic 
times, such as the time to reach maximum flux, and the time scale for the 
decay of streaming anisotropies. Diffusion coefficients have customarily 
been computed from the quasilinear kinetic theory of pitch-angle scattering 
[v. , e.g . , Rowlands , Shapiro , and Shevchenko , I966, Kennel and fingelniann , I966; 

Jokipii , 1966, 1967, 1968, and references contained in Jokipii , 1971] 
along with knowledge of the power spectrum of the electromagnetic 
turbulence. When the turbulence is superimposed on a strong background 
magnetic field, as is the case in the interplanetary and interstellar 
media, the qmsilinear theory contains some well-known difficulties if 
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the Lannor radius of the particles in less than the correlation length 
of the field [ Klimas and Sandri, 1973; Jones , Birmingham and Kaiser , 1973; 
Kaiser , Jones , and Birmingham , 1973; and V<3lk , 1973]« In particular, if 
the magnitude of the magnetic field is constant to first order in the 
strength of the fluctuating field, then peurticle backscattering is not 
allowed. [For a recent discussion of this point, v. Klimas et al. , 1976a,] 
Alternatively, when fluctuations in magnitude are present, the pitch-angle 
scattering coefficient, D , is nearly zero near u = 0 (where is the cosine 
of the particle's pitch angle with respect to the mean field) but at = 0 
is either infinite, in the case of magnetostatic turbulence [Fisk et al. , 
1974 and Gtoldstein , Klimas and Sandri , 1975] ; or, in the general electro- 
magnetic situation, is highly peaked [ Lee and Volk , 1975]. Formally one 
can derive expressions for in this situation [ Jokipii , 1971], but 
interpreting the result in terms of a mean- free path for backscattering, 
which is implicitly done when discussing solar flare events, is suspect 
when scattering through ^i, = 0 is strongly inhibited [ Gal 'per in , Toptygin , 
and Fradkin , 1971 ]. 

Many of the difficulties encountered in the quasilinear theory can 
be avoided if one modifies the trajectory followed by particles to 
include perturbations caused by the fluctuating fields. These perturbing 
forces become especially important for p, ~ 0 because the duration of the 
wave-particle interaction is increased in that region of phase space. In 
recent years several nonlinear theories have been proposed that predict 
significant scattering through p, = 0 [ Jones et al. , 1973 ; V^lk , 1973 ; 
Goldstein, I976] . [An alternative approach that also finds significant 


scattering through m- = 0, while retaining unperturbed trajectories has 
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been developed by Klimas ^ 1976a, b.] In this paper we present 

a detailed comparison between the magnltxide of comxnited from the perturbed 
orbit theory of Goldstein [1976] (Paper I) and that resulting from quasir. 
linear theory. The discussion is restricted to the relatively simple 
situation of magnetostatic turbulence which is a function only of position 
along the mean field (the "slab" model). 

Use of the perturbed trajectory can produce nearly isotropic 
scattering through p, = 0. Then the time for backscattering through 
AO s A (flar ^p.) TT is often about equal to the time for scattering 
through AO ~ 1 and it becomes meaningful to interpret in terms of an 
approach of a particle distribution to isotropy. The value of derived 
from perturbed orbit theory tends to be slightly less than that computed 
from quasilinear orbits. In addition, for power spectra with spectral 
index q ^ 2, the perturbed orbit theory removes the divergence in 

that characterizes the quasilinear result when the distribution function 
is expanded in a perturbation expansion [v. , e.g . , Jokipii , 1966 and 
Hasselmann and Wibberenz , 1968 and I970] . 

In §II we generalize the results of Paper I to include slab models 

with arbitrary spectral indices. We follow the basic approximation of 

that analysis in that we assume, to first order, that the turbulence is 

not modified by the particle distribution (i .e . , a test i>article model). 

In addition, the theory is restricted to the weak coupling limit in which 

it is assumed that tie amplitude of the fluctuating field is sul’ficiently 

small so that interactions between various wave modes are negligible 

compared to wave- particle effects and that the particle orbits are not 
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grossly perturbed during a scattering time— the time necessary to 
propagate %ross several correlation lengths. In $III we evalmte K 

II 

for a variety of parameter values that are characteristic of the inter- 
planetary medlim and discuss the consequences of the results. 
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fin PITCH- ANGLE SCATTERING - A REVIEW 


Ben- Israel et [1975] have recently'- developed a nonlinear 
kinetic theory of strong electromagnetic turbulence. In Paper I their 
work was adapted and modified to describe pitch-angle scattering in 
magnetostatic turbulence in a "slab" geometry. In order to provide a 
comparison with numerical experiments (Kaiser . 1975) » the previous analysis 
was confined to exponential correlation functions to describe the 

statistical properties of the tmrbulence. In such a model, the dlmension- 
leris power spectrum is given by 




( 1 ) 


where, as in paper I, e = X /r is the ratio of the correlation length 

c g 

of the turbulence to the Larmor radius. (Our notation follows Paper I, 
wherein all lengths are measured in Larmor radii and time in Laxmor 
periods . ) 

In this section we first generalize the results of Paper I to include 
power spectra with arbitrary power law indices, q. This will enable us 
to conpute spatial diffusion coefficients from power spectra with 
7/5 i q ^ 2, which is typical of interplanetary observations. For 
mathematical simplicity we generalize equation (1) to the following: 

ifz ‘ 



where f(x) is the gamma function and 2v + 1 = q. The normalization is 
chosen so that 


77: ' 

•km 


( 3 ) 
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where Tr R(r) denotes the trace of the correlation tensor of the field, 
defined by (I-23).* In the slab model, is iiq?licitly defined by 


where RQt) is the Fourier integral transfom of R(r) j and the notation 
follows ( 1 - 3 ^) snd (1-45). Prom (2) eind (3) one has for the correlation 
tensor satisfying (2) fErd^lyi . 1954] 


P(rS • fg-, N 


where X is a dimensionless measvire of length along the mean field 
direction, and where K^(x) is the modified Bessel function of the third kind. 

The set of approximate equations that describe pitch-angle scattering 
in a slab model were derived in §IV of Paper I. From that analysis and 
equation (2) one immediately has for the pitch-angle diffusion coefficient 



where 


<!'?■*'>) casA(h 


* Equation 23 of Paper I . 
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^ ^ rO>*fi) cos Z- (dff cas[jci4>t ’- fiKpf ] 
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s ^^7 Jd\ (Z--X)^ -i- ('^') yC^)^ 


(9) 


( 10 ) 



and T] is the ratio of the root mean square' of the fluctuation magnetic 
field, to the mean field, B^; i«e . , T] a < B'. B'>i/ lB„l. 

Tor convenience we have dropped the subscript "n" on k and Ax. 

For V / ^ the integrals over k in (8) emd (9) must be done 
nuc'^rically, but otherwise the system of equations can be solved exactly 
as before (cf. Appendix E of Paper I). 

For coii5)eLrison, the quasilinear pitch-angle diffusion coefficient 
is given by f Jokipii , I97I; Goldstein et s^. , 1975] 



>j( p(„i ^ 


( 12 ) 


[NOTE : The definition of used here and in Paper I differs slightly 

from that used in Goldstein et al. , 1975.] 
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In Figure 1 we have plotted eolutione of equation (6 ) for 
c ■ 2y 11 ■ 0.1» and v ■ l/3« The results are qualitatively similar to 
those previously presented in Paper I. As before, substantial scattering 
through pi » 0 is present. 


5 m SPATIAL DIFFUSION 


Recently, Earl [1974] has demonstrated t the correct expression 
for X is ;also see Hasselmann and Wibberenz, 1970) 


fri'i 


•I 


( 13 ) 


where s2p,. 

In the past, X has also been computed from the expression [Jokipii, 
II 

1966] 




which should be a good approximation to (13) only for nearly isotropic 
scattering. 

For q ^ 2, equation (13) diverges while (l4) remains finite. 

Earl [ 1974 ] has pre\lously noted that generalizations of quasilinear 
theory which allow scattering through u = 0, also remove the divergences 
in expression (13). 

In the remainder of this section v'e discuss the results of evaluating 
( 13 ) (l4) for various values of T], e, v; and for both quasilinear and 

perturbed orbits. The general conclusion is that use of perturbed orbits 
(equations 6-11) results in approxiinately equal magnitudes of evaluated 
from either ( 13 O or (14). This remains true even for q = 2. V/e dis- 
tinguish the four possible approaches for conputing i^ing a super- 

script QL or ITL to denote quasilinear or nonlinear orbits, respectively; 
and a subscript P or L (perturbation of Legendre) to denote evaluation 

of using (13) or 'l4), resi>ectively. 
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We confine our attention to values of c between one and three euid 
T) ^ 0«6« At small c (high energies) the slab model becomes ruiphysical, 
but in more realistic turbulence geometries the usual quasilinear 
treatment is adequate because the unperturbed trajectories are straight 
lines and consequently the difficulty for l^rge c of treating ^ = 0 
disappears f Jokipil , I966, I967, I968; Klimas and Sandri , 1971]. At low 
energies (c » 1) problems arise because in the interplanetary medium T) can be 
0.6 or larger. In this case equations (6) - (11) become inaccurate 
represen'' ations of the perturbed orbits for c > 3 and *». complete 
derl'vatlor of the perturbed orbit equations has not been carried out 
for arbi^4•az•; *</alues of u* [See Paper I for a more de'bailed e^valuation 
of nonlinear effects when p, = 0.] In addition, we show in the Appendix 
that the weak -coupling approximation, which is essential to the deri-vation 
of (6) - (11), breaks '’own for T) c ^ 1. 

The dimensionless results can be easily transformed into more 

familiar dimensional units as follows : Let X denote the dimensioncl 

II 

diffusion coefficient, then 


when) ud B.* <^8^ 

which we take to be 5 x 10 ^ G in the interplanetary medium. 


In Tables I and II we show K for T| = 0.1, e - 1, 2 and 3; and 

II 

V = ^ and 1/3 (q ® 2 and 5/3). Use of perturbed orbits has completely 

remo-ved the dl 'mergence in for q = 2 .vnd yields 

II II ll»I' II 

For q = 5/3 t-»e four forms of K are all nearly equal. The success of 
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the perturbed orbit theory in removing the divergences present in 

quasilinear theory for steep power spectra (q s 2) is, of course, 

common to any approach that gives a non-zero value to D at = 0. 

IJ' 

The near equality of to K^_ arises becaiise scattering through 

II II 

[i 0 is nearly isotropic, and eith:*r spnerical harmonic or Legendre 

expansions of the distribution are expected to yield good approximations 

to K . 

11 

In the interplanetary medium T] is rarely as small as 0,1. However, 

it can be determined from power spectra computed from data collected 

from magnetometer experiments on space probes. We discuss below 

two examples of such spectra; the first from Jokipii and Coleman [I968] 

(Figure 2) and the second from Fisk and Sari [1973] (Figure 3) . From 

Figure 2, v =** 0.2 (q 7/5) while from Figure 3j ~ l/3 (l ^ 5/3) and 

we assume, for illustrative purposes, that the spectra shown represent 

a slab model. P(f) is the power spectral densit'-- in y jYa. (1 y = 10 ''^g) 

as a function of f = V k/2TT, where V is the speed of the solar wind, 

sw sw 

taken to be 350 km/s. P(f) is related ^(k||) through the definition of 
R(X) ; viz . 

iCF = 2 cos A-J 


(16) 


SO that 






cff -p(-n 


(17) 
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where 


P(f3 H V(^X.VVs^ 


Equation (17) also defines the total variance of the fluctuating field; 

2 

which, for Figure 2 is 4.3y (T1 = .6), (the Fisk and Sari spectrum is 

11 

discussed below). For the spectrum in Figure 2, 2 x lo cm. The 

resulting values of are shown in Table III, and again use of perturbed 

orbits does not significantly modify . 

Fisk and Sari [1973] have argued that low-energy particles 

(<^ 100 MeV/nucleon) do not interact with tangential discontinuities 

that are included in computing the power spectrum shown in Figure 2. 

The authors th t go on to remove tangential discontinuities from their 

data set and compute the spectra labeled in Figure 2. In this 

10 

case T] a: 0.3 (Sari, private communication), and ^ 2 x lO"" cm. The 
resulting values of K are shown in Table IV. 


The values of K in Table IV (and to some extent, also those in 
II 

Table III) are consistently smaller than those deduced from the decay 
portion of solar flare events [v. , e,g . , Webb ^ al., 1973; and Countee 
and Lanzerotti, 1976] . Counbee and Tanzerotti [I976] conclude that an 
empirical fit to the particle data gives 2 x iO cm /s (at proton 

energies of ~ 7 MeV) , a factor of ~ 7 larger than shown in Table IV. 

Webb et al. [1973] also found that the scattering mean free path for 
various solar particle events was smaller than could be understood 
theoretically from numerical solutions of the Fokker- Planck equation 
which describes spatial transport. The resolution of this discrepancy i 
really not well under s tooc' . -ii tinere are several possibilities, some 



of which we mention below: First, we have restricted the turbulence 


model to one in which all the k- vectors are parallel to which tends 

to underestimate K . Volk [1973] has argued that the radial 

II 

expansion of the solar wind causes the wave vectors to refract into the 

radial direction, thus increasing K . Recently, Morfill [1975] and 

Morfill et al . [1976] have investigated this possibility in more detail, 

with the conclusion that one can indeed ^ increase K by a much as a 

II 

factor of 10 if the wave vectors of interplanetary Alfv6n waves are 
predominately in the radial direction. Their analysis includes a dis- 
cussion of effects due to an admixture of compressive modes. Unfortunately, 
there is little definitive information about the direction of propagation 
of interplanetary magnetohydrodynsimic (MHD) waves and no evidence has been 
presented that the wave k- vectors are in the radial direction. Second, in 
order to fully understand interplanetary propagation it is undoubtedly 
necessary to have a fairly detailed knowledge of the nature of the 
magnetic turbulence before, during and after a solar flare event. Such 
a detailed study is not generally made, and consequently is not often 
evaluated for the actual turbulence at the time of a flare. There are, 
however, several analyses which have related to the magnetic power 
spectrum at the time of the particle events (v., e.g . , Lanzerotti et al. , 
1973; VJebb ^ £l., 1973; and Wibberenz et ail., 1973? 1976). A general 

conclusion of these studies is that K evaluated from the weak-coupling 

II 

theory is systematically smaller than implied by the observed time 

behavior of the particle fluxes. Third, one does not know the heliocentric 

dependence of K . If (Sari, 1973) ^ increases rapidly toward the Sun, 

II II 

then large mean free paths inside of 1 A.U. could explain the flare 

lU 


observations at 1 A.U.; however, large values of are also needed to 

understand the small gradient in cosmic ray intensity that is observed 
¥ 

on Pioneers 10 and 11, so again one cannot completely remove the dis- 
crepancies in this way. One could also argue that if ““ 2 x lO^cm, 
then 10 MeV protons correspond to e 20. Then, with T| “• 0.6 one has 
seriously violated the weak-coupling approximation; and, as there is no 
extant kinetic theory valid for strong coupling in strong turbulence, 
one might expect th*»t a more correct theory, when it appears, would 
improve the situation. Even this explanation is seriously weakened, 
however, if Fisk and Sari [1973] are correct that low-energy particles resi)ond 

'’between '’eEAL- '‘"d e “ 2 

(12 MeV protons) the weak-coupling approximation is not badly violated; 

but from Table IV, one still finds K to be too small. Yet another 

II 

approximation that is always made in deriving is that the particle 
distribution behaves as a collection of test particles. While this is 
an excellent approximation at moderate energies it may well breakdown 
at the low energies (< 10 MeV for protons) discussed by Countee and 
Lanzerotti [1976]. A crude estimate of the energy density of the particles 
indicates that at times it equals the energy density in the fluctuation 
fields; a situation which in principle necessitates a full self-consistent 
treatment of the Vlasov-Maxwell equations for a complete theoretical 
treatment. [Such a theoretical framework is available only in the 
weak-coupling regime (Ben- Israel et al. , 1975)*] Finally, Nolte and 
Roelof [1975] have developed a mathematical treatment which examines the 
consequence of assuming ’’scatter- free" propagation at the onset. 

iirppwODUCffiiur/ of the 

PAGE IS POOR 
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5V CONCLUSION 


We have computed the spatial diffusion coefficient using both 

perturbed and quasilinear orbits within the context of the weak-coupling 

approximation. Use of the perturbed trajectory removes divergences 

found in quasilinear calcxilations for steep power spectral indices 

(q i 2) when equation (13) is used to conpute . For less steep spectra, 

l<a<2,K«^ tends to be slightly less than K ^ due to enhanced pitch- 
11 II 

angle diffusion through ij, = 0, If the correlation length of the inter- 
planetary turbulence is 2 x 10^ cm, then for e < 3 (proton energies 
> 4CX) Mev) one obtains values generally consistent with observations. 

At lower energies (e » 3), the weak-coupling approximation is violated 
and one may require a kinetic theory valid for strong co\pling. 

However, if Fisk and Sari [1973] are correct in their assertion 
that low-energy particles (*“ 30 Mev protons) traverse a stochastic 
field characterized by 2 x lO^^cm, then one should, in principle, 
be able to use the present theoretical formalism to determine the spatial 
diffusion coefficient for ~ 30 Mev protons. The resiilts of such a 

computation, using either perturbed or quasilinear orbits, yields spatial 
diffusion coefficients that are consistently too small by factors of 5-10. 
[Unless there is significant focusing of the wave k- vectors into the radial 
direction.] However, it is possible that for this low-energy particle 
population the test particle approximation has broken down. Correction 
of this breakdown would require significant modification to the theoretical 
formation that is now used, even within the context of the weak-coupling 
approxima.tion . 
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FOR 0,h TO 2 GeV PROTONS (q = 2 and T] = 0.1) 



FOR 0.4 TO 2 GeV PROTONS (q = 5/3 and T] = 0.1) 
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FIGURE CAPTIONS 


Figure 1. 


Flg\ire 2. 


Figure 3. 


The dimenslonlesfl pltch>ongle diffusion coefficient, e 
is plotted as a function of The solid line is the 
result of using eqmtions (6) - (11) (nonlinear orbits) , 
while the dashed line results from use of equation (12), 
enjoying unperturbed orbits (quasilinear approximation). 
The spectral index is q = ^/3, and 11 and e are 0.1 and 
2,0, respectively. 


Power spectrum of the component of the interplanetary 
magnetic field normal to the solar eqmtional plane, 
observed on Mariner 4 . A fit to these data with an 
analytic function of the form of equation (2) yields T) “■ 0.6 
and 2 X 10^ cm (after Joklpii and Coleman , I968) . 


Representative spectra for the total observed power density 
(Pj^l)» power density due to directional discontinuities 
(Pjjlg), and the powe^ density due to field fluctuations 
between the discontinuities (^3g«jy2EN^* ^ these 

data with an analytic function of the form of eqmtion (2) 


yields, for 


0.3 and X «»• 2 X 10^° cm (after 
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APPENDIX 


The weak- coupling approximation, as used here, is essentially the 
claim that (Appendix B of Ben- Israel et al, [ 1975 ] or 1-13) 




(A-l) 


This is true whenever 

.r 




(A-2) 


All quantities in the integrand are 0(1), and the integral converges in 
after a time T^lorT^e(i.e., after a cyclotron period or a 
correlation time, which are the two characteristic t: ae scales of the 
problem. Therefore, (A-2) is true if T] « 1 and T] e « 1. (This is 
identical to the conditions of Ben- Israel et al.[i 975 ] that CXg « 1 and 
ctg « 1.) For low energies (e » 1) one expects (A-2) to saturate 
fter many Larmor periods (x > e) and therefore one needs T] e « 1 for 
(A-l) to hold. 

Similarly, the perturbed trajectories are likely to be a good 
approximation :i f (cf. equation B3 of Ben- Israel et al . [I 975 ] or I-I 6 . 

uctr)= ULTr) 

with 

n jM \ ^'ua) > ] <^ / I'-'O 


Iteration leads to 

Ji aitr~\ ) jf'ucky « ( 


which is true if T] « 1 and T] e « 1 j as above. 
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